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ABSTRACT 



The general equations of motion for a rigid body are 
derived in cylindrical coordinates by Lagrangian dynamics 
and used to model the motion of an airplane in a steady spin. 
After simplification, the equations are cast into a form 
utilizing conventional aerodynamic derivatives along with 
other derivatives which may be significant in spins. An 
iterative numerical solution procedure is outlined which 
should simplify the problem of solving the nonlinear dif- 
ferential equations, and relationships between the Euler 
Angles used in the equations and the more familiar ordered 
set of pitch, roll, and yaw are derived to permit computer 
input and output of orientation to be more easily visualized. 
The inverse problem is also considered, and equations are 
derived to translate spin test data into the six Lagrange 
coordinates used so that computer calculations may be checked 
against test results and actual aerodynamic forces can be 
compared with computed values. The paper concludes with a 
discussion of the implications of this spin analysis along 
with applications and extensions. 
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TABLE OF SYMBOLS 



The definitions 


of the symbols used in this paper are as 



follows ; 



C9 


Center of gravity 


c 


Mean aerodynamic chord, feet 




Aerodynamic force coefficient along the 
subscripted axis , dimensionless 


Cig 


Aerodynamic moment coefficient about the 
subscripted axis, dimensionless 




Aerodynamic force along the subscripted 
axis, lbs. 




Lagrange generalized non-conservative 
force, lbs. 




Lagrange generalized force in the subscripted 
coordinates, lb. or ft - lbs. 


9 


2 

Acceleration of gravity, ft/sec 


9jf , 9y-» Si 


Airplane eg acceleration components along 
the subscripted axis, ft/sec2 




Airplane principal moments of inertia, 
slug - ft^ 


M 


Airplane mass , slugs 




Moments about the subscripted axis, ft - lbs. 


»? 


Sequential subscript, integer 




Angular rates of rotation about the airplane i 

X, Y, and Z principal axes, respectively 




Lagrange generalized coordinate 




Spin radius, feet 


6 


2 

Airplane wing area, ft 


T 


Lagrange kinetic energy term, ft - lbs. 


t 


Time, seconds 
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At 


Increment of time, seconds 


V 


Lagrange potential energy term, ft - lbs. 


% 


Airplane eg velocity ft/sec 

Airplane eg relative wind velocity component 
along the subscripted axis, ft/sec 


y. ? 


Airplane principal axis system 




Reference Cartesian axis system 


2. 


Airplane altitude, feet 


« 

^0 


Rate of ascent, ft/sec 




Airplane angle of attack relative to principal 
axis system, degrees 




Direction cosines in terms of Euler angles 
(see Table I) 


/^ 


Airplane yaw angle, degrees 


y 


Lagrange generalized coordinate indicating 
direction of spin radius vector, radians 


y 


Spin rate, radians/sec 


^4.. ^»-y 


Aileron, rudder and elevator deflections, 
degrees 


(t> 


Airplane principal axis system orientation 
Euler angles with respect to the Xj^, , 

system, radians 




Airplane orientation in pitch, roll, and 
yaw with respect to the system, 

radians 


®P> ^y 


, Pitch, roll and yaw rates, radians/sec 


^xi/ y ^i( 


Incremental angular rotations about the X^^ , 
^1' ^1 system axes, radians 


f 


3 

Air density, slugs/ft 




Angular rates of rotation about the air- 
plane principal axes, equivalent to p, q, 
and r, radians/sec 
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I. 



INTRODUCTION 



A. FOREWORD AND ACKNOWLEDGEMENTS 

The intent of this work was to develop an alternative to 
the conventional method of computing spin characteristics 
which thus far has not been uniformly successful in pre- 
dicting airplane spin modes. In this effort, the author 
has tried to approach the problem with the capabilities of 
modern computers in mind and devise a method of solving 
directly for spin modes rather than mathematically tracking 
airplane motion through all intervening phases of motion 
preceding the spin. If any success is realized through this 
approach as outlined herein, it will be due in no small 
measure to the tireless efforts of Lieutenant R. L. Champoux, 
USN, whose ccatribution of months of hard work and extensive 
knowledge of computer programming were instrumental in the 
evolution of this approach. Lieutenant Champoux 's separate 
work on the computer aspects of this analysis may be found 
in Ref. 9. The author is also indebted to Professors L. V. 
Schmidt, T. H. Gawain, and M. H. Redlin, of the U. S. Naval 
Postgraduate School whose counsel helped shape many aspects 
of this work. 

B. BACKGROUND 

Airplane spin or autogyration is a mode of motion 
exhibited by virtually all fixed wing aircraft. It is at 
best an undesirable flight maneuver and, with growing 
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frequency in new airplane designs, a terminal maneuver from 
which recovery cannot be effected. Although the flight con- 
ditions which lead to spins can be avoided, pilots of military 
light attack airplanes and especially fighter airplanes do 
not have that option if pressed into flying to the "edge of 
the airplane's flight envelope" by the mission or a tactical 
situation. 

The need to utilize the edge of the flight envelope in a 
combat role is readily apparent when one considers that it is 
only in this fringe area that the superiority of one airplane- 
pilot combination over another can be realized. However, less 
obvious is the extent of additional exposure necessary for 
training in this area if a pilot is to become proficient in 
a particular airplane under such conditions. Thus, flight 
near the edge of the envelope must be more than an occasional 
event for the military pilot. Unfortunately, loss rates due 
to irrecoverable spin modes in advanced tactical jet airplanes 
frequently makes effective training in this area economically 
prohibitive. Thus the margin of superiority built into the 
airplane at great expense is in many cases lost. 

It is in recognition of this fact, and the fact that spin 
predictions based on conventional analysis techniques have 
not shown much success, that this effort was undertaken. It 
is not expected that this analysis of airplane spin will lead 
directly to solutions of the spin problem itself, but it is 
hoped that this work will cast the problem in a new mathe- 
matical perspective, from which additional insight into the 
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spin problem might be gained towards a goal of eventually 
designing out irrecoverable spin modes. 

C. GENERAL 

Historically airplane spin studies have employed forms of 
the same equations of motion which are used in the small 
perturbation approach to airplane dynamics. This method of 
mathematically modeling the spin appeared promising because 
the equations could be readily integrated by numerical 
methods to produce time-histories of characteristic parameters 
which could then be compared directly with flight test time- 
histories. It was also apparent that the body axis coordinate 
system used in the classical approach provided simpler 
equations for random tumbling motion than those obtained 
using an inertial system of coordinates. Thus the post- 
stall and incipient spin gyrations of the airplane, which in 
many cases are of equal concern, can be more readily tracked 
by a computer in a body axis system than an inertial system. 
The degree of success realized by thus extending the classi- 
cal approach however depends to a great extent on the 
judgment and intuition of the investigator in leaving just 
enough coupling in the equations to permit satisfactory 
representation of the motion without introducing excessive 
non-linearities. Frequently this process degenerated to a 
trial and error parametric approach without much insight 
being gained and without general validity being obtained for 
all airplane configurations or spin modes. 
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On the other hand , the treatment of the spin problem in 
an inertial coordinate system appears to add undue complexity 
to the computer tracking problem, and indeed for random 
tumbling motion the non-linearities become excessive. However, 
if the steady, or quasi-steady^ spin is to be studied, an 
inertial cylindrical system offers a number of distinct 
advantages, among which are the ability to differentiate 

2 

between motion due to spinning and motion due to tumbling , 
and a sufficient reduction in the complexity of the coupling 

3 

to permit computer solution of the exact equations . Since 
the intent of this study was to develop a generalized method 
of solving for steady spin modes, the equations of motion 
were derived in an inertial coordinate system. The details 
of this coordinate system are shown in Figure 1, and discussed 
in some detail in Section III. 



Spinning motion in which state variables have a mean com- 
ponent and an oscillatory component. 

2 

Random motion inconsistent with steady or quasi-steady spin 
motion. 

3 

No small quantity approximations are made in either the 
derivation or simplification of the equations of motion. 
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II. SCOPE 



This work was preceded by a survey of current airplane 
spin literature [Ref. 1] to determine the present state-of- 
the-art in airplane spin analysis. From that study it was 
evident that possibly another avenue of approach to the 
problem ought to be explored. Thus, rather than attacking 
the spin problem on the broad formidable front of mathe- 
matically following an airplane through all transient 
gyrations from stall to spin, it was decided to limit the 
scope of the investigation and probe directly into the final 
steady state motion of the spin where a force and moment 
equilibrium situation might yield to simpler mathematical 
representation. In this approach no restrictions were 
imposed on the form of the airplane aerodynamic data to be 
used in this analysis other than requiring that it be ref- 
erenced to the airplane's center of gravity and principal 
axes. By leaving the details of the makeup of such data 
somewhat general, the analysis is freed from any incumber- 
ances of classical aerodynamic derivatives and thus can 
proceed more directly to the derivation of the general 
equations of motion. Although conventional force and moment 
derivatives are suggested for use in initial computer pro- 
gramming, subsequent work may suggest a more appropriate form 
of aerodynamic force and moment data, or possibly more 
representative aerodynamic derivatives may be devised. 
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III. DERIVATION OF THE EQUATIONS OF MOTION 



The coordinate system used in the derivation of the 
equations of motion is a combination of an inertial 
cylindrical system and the airplane principal axis system. 

The vertical (z) axis of the cylindrical system corresponds 
to the central axis of the motion or "spin axis". The 
advantage of the cylindrical system is that it locates the 
airplane center of gravity in convenient terms of an alti- 
tude coordiante (z^) / spin radius coordinate (R) and angular 
position coordinate (?) , when the airplane is in a steady- 
state spin. The orientation of the airplane about its 
center of gravity with respect to the inertial cylindrical 
system is specified in terms of Euler angles. A cartesian 
coordinate system is fixed at the eg position on the (R) 
vector having its (x^^) axis in the (-R) direction and its 
(z^) axis parallel to and in the opposite direction of (z^) . 
This cartesian system provides a reference for the orientation 
Euler angles. Thus for zero values of the Euler Angles (0), 

(y) , and ((p) , the airplane would be upright, wings level, 
and have the (X) axis of its principal axis system pointed 
inward, along the (R) vector, directly at the spin axis. 

The details of this hybrid coordinate system are shown in 
Pig. 1 along with the positive directions of the position 
coordinates. A general schematic of Euler angle relation- 
ships is contained in Fig. 2 , and equations for the 

direction cosines of the airplane axes expressed in terms of 
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J 





Fig. 1. The Coordinate System Used to Model the Spin 




Fig. 2. A General Schematic of Euler Angle Relationships 
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Euler angles are listed in table I. The rationale behind 
this choice of coordinate systems is that it more simply 
(in a mathematical sense) represents the motion being 
modeled. Thus even for oscillatory spins the motion would 
be readily recognizable and translatable into mean values 
of the six coordinates (three position and three orientation) . 
In addition, since the principal objective is to solve for 
steady spin modes, the complexity of stall and departure 
motion in this coordinate system has no impact on the results. 

4 

Table I: Direction Cosines in Terms of Euler Angles 



Cosines of Angles between X,Y,Z and Xl,Yl,Z^ 


• 


X 


y 


Z 


Xx 


flu = cos ^ cos ^ 

— sin 0 sin cos e 


fl2i = — sin <f) cos ^ 

— cos ^ sin ^ cos $ 


fl3i = sin e sin 0 


Yx 


®X 2 “ cos <p sin ^ 

4- sin cos ^ cos $ 


®22 “ sin <f> sin ^ 

+ cos 0 cos 0 cos e 


'X. 

fl32 = ' — sin $ cos 0 


Zx 


cti3 = sin ^ sin ^ 


fl23 ~ sin $ cos 0 


fl33 = cos 8 



With the coordinate system defined, the equations of 
motion may be expressed in terms of the six generalized 
coordinates z^, R, , 0, ^ and <f> by the use of Lagrangian 
mechanics. The only simplifying assumption made in the 
derivation of the equations is that the airplane develops 
zero thrust in the spin. This is frequently the case in an 
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Reference 2, page 158. 



actual spin situation due to jet engine stalls and/or flame- 
outs resulting from inlet airflow distortion. However, the 
restriction imposed by this assumption can be easily removed 
by including negative dissipative work terms on the right 
hand sides of the Lagrange equations which follow. Such 
terms would include X, Y and Z force component's as well as 

force and moment derivatives about all three axes due to o( 

* • 

and ^ ^ and jS , and the angular rates p, q, and r. 

The derivation of the equations of motion is best illus- 
trated by first writing the general Lagrange equation: 




Where : 



T = System kinetic energy 

V “ System potential energy 

q^ “A generalized coordinate, r = l,2,**-,6 

Fq^ “ A generalized force due to non-conservative 
forces, r = 1,2, ••*,6 

For this particular coordinate system, T can be written as: 



Where Ig; / * 60y » CO 2 ' moments of inertia and 

angular rates of rotation about the airplane's principal 
axes. 

Using Figs. 1 and 2, the angular rates relative to the 
inertial reference can be written in terms of the spin 
coordinates as follows : 
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'i-ecos(p 



(3) 



-esitjp 



(4) 



co^ = )o(^2> "^4^ 



(5) 



Note that f is included to account for the rotation of the 
Y^, system. 

The partial derivatives for 

and (J) are : 



SJ„ 



( 6 ) 



« M 

"Z * 

57 ^^^3 '^y6j/V23 ■*■ <^^3 

^ = XyO^CoSp -lyCOySi^Jp 

* Xx + ly 63/ «^Z3 + X^ Oj «<33 



(7) 



( 8 ) 



( 9 ) 



( 10 ) 




1 ^ 



( 11 ) 
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I 

J: 





And the time derivatives of these are ; 

•^Iy(C0y<<iS+^y<^23) ^i(^g <?i ''' <^e ) 



Ji ^ Ij) "" (^i^CoS (p-(pQx StKj <j?)-Jy((OyS/^^-h<l>(A^yC(/S(l>) 
+l5 (O-^i 0(33 ‘5<'jj) 



Next, the partial derivatives are obtained 









I? - 






(12) 

(13) 

(14) 

(15) 

(16) 

(17) 

(18) 

(19) 

( 20 ) 
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Se 



5^ 



g=I.CO.‘^^-.W^.X,Ci^^ 

And from V “ MgZ the values of are 

o 'Qfr 



( 21 ) 



(22) 



(23) 



(24) 



~ ^ - *C>1/ „*5l/ . (25) 

Next, evaluating the quantities » ^/3 

where k “ x, y or z and i ~ 1, 2, or 3 , and r = 4, 5, or 6 : 



-c^z. + eccs<P -eh‘'^4> 



(26) 



= C^V y ,) 0 ( 2 ^ <k2i-9StNp^eif> cos p 



(27) 






(28) 



0 

= e>cosesiN<j) e cos<p 



(29) 



ci'2S-GC0SOC0S<P - <p ^ S/AJ<p 



(30) 
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0^33 = -esiKie 



(31) 



'do 



^3 = Cose SiKXp 



(32) 






^ S/a; 0 CoS 



- Cos 0 Cos <p 



' ^0('23 

^<p 



^ - siKi e Sffj <p 



(33) 



(34) 



(35) 



7 e" 



^ - 5ito e 



(36) 



'3o(', 3 *^^23_^^33 po/33 ^ 

The partial derivatives , etc., are also easily obtained, 

and will not be shown here. 

Finally, by substituting the derivatives defined in 
equations (12) through (25) , six Lagrange equations of motion 
can be written in terms of the quantities defined in equations 
(26) through (37) and the partial derivatives of the angular 
velocities are as follows ; 



A)?, = Fg, 



(38) 
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(39) 



Af/e'-Afr^/e*F^ 

Ix ((^ ^)3 ) + *3^y + CO;. c^r j 

+X 2 C^02«<33+ C02o;r3^j = Fy (40) 

”Xy(c0y«5/^ (p-tcOy^ COS^) 

-(tyU, + Ix ^Ox + 1^«)^ II?; = Fe (41 ) 



Xx(^X®C,3 f- C0> -f 2/ (cby ^2S 4 CO^ C^2J ) 

• +X2 ( oc'^^ + ‘^33^ “ ^ (42) 

= F 54 («) 

As discussed in Ref. 2, the concept of virtual work can 
be employed to derive the non-conservative forces F .fp . 

Ff I Fq , F^ , and F(p which constitute the generalized forces 
in the Lagrange equations (38) through (43) . In this 
technique, each of the generalized coordinates (Z^/ R/ ^ / 

Q , ^ or P ) are increased, one at a time by a small positive 
increment while holding all other coordinates fixed, and the 
expressions for the virtual work terms are determined. Since 
system forces and moments due to aerodynamic effects are 
applied at the eg and expressed in orthogonal X, Y coordinates, 
no coupling of the virtual work terms occurs and they can be 
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readily derived by inspection using pigs* 1 and 2. The 
partial derivative of the sum of the virtual work terms is 
then taken with respect to one of the generalized coordinates 
and the corresponding generalized force is thus obtained. 

For example : 




(44) 



Where , and are positive incremental angular dis- 
placements about the respective axes , and 



Since, by inspection of ig. 1; 




(45) 



And, as a result of orthogonality: 




(46) 



Thus we obtain: 



-Fi 



(47) 



And, by identical procedures: 






(48) 




(49) 



^ c Cos My, ^/AJ ^ cos^ +My S/fJS 



(50) 
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( 51 ) 



p - Mf, SnJe S/*J - My, S/fJ 0 as f 
+ Mi, cos © = 



Now, by employing the direction cosine relationships listed 
in Table I these forces and moments can be obtained in terms 
of forces and moments along and about the airplane's principal 
axes. Thus equations (47) through (52) can be rewritten as; 

^z/ ^ ^3/ ^z) (54) 

- Mx ^^3 - My S/A/(p ( 56 ) 



= ^2 



(57) 



icx:,}Mx +^2S^Y 



(58) 



Finally the complete dynamic equations governing motion in the 
chosen coordinate system can be written in terms of the six 
generalized coordinates z^ , R, , Q,^, and (f? by equating 
the left hand terms of equations (38) through (43) with the 
respective right hand terms of equations (53) through (58). 
These equations are written in their entirety in Appendix A. 
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IV. INCORPORATION OF CONVENTIONAL AERODYNAMIC DATA 



A. DERIVATION OF REFERENCE PARAMETERS 

In order to utilize conventional aerodynamic force and 
moment coefficient data in the equations of motion just 
derived it is necessary to obtain such typical reference 
parameters as airplane angle of attack {<K) , sideslip angle 
if3 ) , center of gravity velocity ( 1/^^ ) , roll rate ( p) , 
pitch rate (^), and yaw rate (T) in terms of equation 
variables and their time derivatives. This can be done with 
the aid of Figs. 1 and 3 as follows: 

Vc^ = +y, Cs (59) 

where unit vectors along the Xj^ f / 

axes respectively. 

But from Fig. 1: 

Thus : 

-R ^5 (60) 

Therefore the relative wind at the airplane eg is: 

^ (61) 

By use of the direction cosines given in Table I, the relative 
wind can be resolved into components along the airplane's 
X, Y and Z axes as follows: 
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(63) 



^ '^^22 y /^ 

VfZiA/g^ = C<3i ^ ’*' ^Jz. 



1^ Y 




tan C>( - 



V 

RWZ 

V 

RWX 




V 

RWY 

V 

RWX 



Fig. 3 Definition of Positive ^ and ^ Angles 
Therefore and can be written as : 




(65) 



^ f o(2ik \ 

The magnitude of the relative wind can be written from 
equation (61) as: 



* ("W^+ z/'j 



(67) 
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Finally, the rotational rates p , q , and r about the air- 
plane's X, y , and z axes are identical to the angular 



B. MODELING THE FORCES FOR THE STEADY SPIN 

Tabulated aerodynamic data, although virtually always 
appearing in the form of dimensionless coefficients, varies 
widely in terms of the reference axis system, point of 
application of forces and moments, and the characteristic 
length used in non-dimensionalization , depending on the 
intended application of the data. By the use of the princi- 
pal axis system in the foregoing derivation of the equations 
of motion, considerable simplication was achieved. There- 
fore it appears that conversion of aerodynamic data to the 
principal axis system and eg as references would generally 
be preferable to the alternative of adapting the equations 
to the particular form of the data. This should present no 
serious problems, however, since the equations must be 
solved numerically and conversion of the data by computer 



rates , COy and CO, given by equations (3) , (4) , and 



(5), and can be written directly as: 



P^( f e cosep 



( 68 ) 



sif^ ^ 



(69) 




(70) 
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is by comparison a relatively minor additional task. Since 
some data manipulation will likely be required in any event, 
the mean aerodynamic chord ( C ) has been used in the follow- 
ing development as a common characteristic length in an 
attempt to be consistent and also allow the use of a single 
"characteristic" length parameter as a computer program 
input. However, the equations which follow could be adjusted 
appropriately to use other characteristic lengths. 

The general form of the aerodynamic forces and monents 
contained in the right hand terms of equations (53) through 
(58) can be expressed as follows: 



Where C has been used as a characteristic length in non- 
dimensionalizing all moment coefficients. The force and 



Fi = '/z <^f=; I ‘ X, r, ='■ i 

Ml = '/xC»;ev‘Sc 



moment coefficients and Cm; ) are, in turn, of the 

following general form; 











Where : 



1. All Cfjj and jj ( j “ 1 , 2 , • • • 7 and i ~ x ,y or z) 
are aerodynamic force and moment derivatives, as 
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defined in Table II below, which must be experi- 
mentally determined and tabulated for a wide range 
of p< and ^ values. 



Table II; Definitions of Aerodynamic Derivatives^ 





SUBSCRIPT j 


i“x ,y , 
or z 


1 


2 


3 


4 


5 


6 


1 


^Pii 








3Cfi 




9Cp; 


dCn 


m) 








dir 














3cmi 






CnVj 


m) 


m) 






dir 





2. -f(P), 'f(i), and f (r) are functions of roll, pitch 

and yaw rate respectively. 

3. ^4 ' defined as : 

~ Aileron deflection 
“ Rudder deflection 
^Ij = Elevation deflection 

A good example of the form of tabular data that could be used 
in this manner is contained in Ref. 3. 

At this point some comments are appropriate to insure 
that the above method of handling the aerodynamic forces 
and monents are put in proper perspective. First of all it 
should be noted that the validity of any calculations of 
spin parameters obtained by this approach to the spin 



Aerodynamic derivatives with respect to and have not 
been included since and ^ are zero for the steady spin. 
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problem depends entirely on the modeling of the aerodynamic 
forces and monents since the equations of motion as derived 
in Section III are exact. Secondly it should be recognized 
that the aerodynamics of the spin problem involve low speed 
three dimensional stalled flow, and as such forces and moments 
are nonlinear. Thus the superposition of the separate 
effects of triaxial rotation, control deflection, and other 
static wind tunnel data cannot be justified analytically. 

The only wind tunnel test technique that presently appears 
to offer any possibility of yielding satisfactory aerodynamic 
data for spin is the rotary balance force and moment measuring 
apparatus which is reported on in Ref. 4. 

Next, allowing that for want of better data the conven- 
tional approach outlined above will be used to exercise 
computer program and look for a zeroth order approximation 
of spin characteristics, there is another problem that stems 
from this application of conventional data. The problem is 
that existing data is generally of insufficient range in o( 
and to accommodate the high angles of attack and side- 
slip angles characteristic of autogyration. Also many of 
the basic set of 42 aerodynamic coefficients suggested by 
the above development describe coupling effects which in 
conventional airplanes at controlled flight altitudes are, 
in general, not significant and therefore have not been 
measured. Thus, for example, roll moment due to elevator 
deflection is normally ignored, but may be of some signifi- 
cance in a spinning F4 airplane where the strong negative 
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dihedral of the stabalator could cause the upwind panel to 
act as a rudder with prospin deflection at highc^ and^ 
angles. Certainly in such cases none of the aerodynamic 
derivatives including the 42 above can be ignored without 
testing. 

Now admitting to the fact that force modeling with con- 
ventional stability derivatives will be done with at best 
an incomplete set of data, it then becomes necessary to 
exercise caution to insure that those coefficients that are 
available and used will allow a force and moment balance to 
be somehow achieved. Thus, if at least one damping term is 
not included for every driving term in each of the six 
degrees of freedom, obviously no stabilized solution is 
attainable. Finally, in fairness to this effort it must be 
pointed out that the modeling of aerodynamic forces and 
moments on the spinning airplane, although recognized as 
the key to the problem under study, was non-the-less 
peripheral to the expressed purpose of this preliminary 
work. Once a workable computer program has been written 
and checked by use of conventional aerodynamic data, sub- 
sequent efforts should be directed towards improving the 
modeling of aerodynamic forces and moments. 



31 



V. SOLUTION OF THE EQUATIONS OF MOTION 



A. GENERAL CONCEPT 

From an inspection of the equations of motion contained 
in Appendix A (in abbreviated format) it is evident that, 
although they are ordinary second order differential 

equations, they are also so extensively coupled and non- 

\ 

linear that direct analytical integration is out of the 
question. Even if the left hand sides of the equations 
were linear, the forces and moments which constitute the 
right hand sides are functions of a three dimensional flow 
involving the geometry of the airplane. Consequently any 
solution scheme will of necessity involve some form of 
iterative numerical procedure to circumvent direct inte- 
gration. The problem thus becomes one of finding an 
iterative scheme which provides some degree of uniform con- 
vergence to steady or quasi-steady values of the six spin 
parameters { ^ , y , 0 , ^ , and p ) without losing 

numerical significance or using excessive computer time. 

If the problem is visualized as a trial and error 
positioning and orientation problem where the airplane is 
incrementally reoriented angularly in response to an aero- 
dynamic moment inbalance and simultaneously adjusted in 2^, 
/and in response to an aerodynamic force inbalance, it 
then appears to be reducible to a simple series of stepwise 
adjustments or iterations. Also, since the problem would 
thus be restarted on each successive iteration with updated 
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force and moment data corresponding to the adjusted spin 
parameters, such a procedure would avoid the possibility of 
loosing numerical significance. The problem of determining 
consistent incremental changes in each spin parameter to 
achieve some degree of uniform convergence of such an 
iterative scheme can be easily handled by using the principle 
of conservation of momentum. By rearranging the equations 
of motion in such a way as to facilitate the computation of 
accelerations in each of the six spin coordinates it is then 
possible to apply a common increment of time and thereby 
compute a dynamically consistent set of incremental changes 
in the parameters and obtain new parameter values for the 
subsequent iteration. This, in a sense, amounts to follow- 
ing the motion of the airplane as is done in the conventional 
moving coordinate system approach to the problem, but the 
hybrid inertial coordinate system has several distinct 
advantages in this application, as discussed below. 

The rationale upon which this suggested solution techni- 
que is based is that it duplicates in a computational sense, 
what actually happens to an airplane under such flight con- 
ditions. However one key advantage that the computational 
scheme has over the actual airplane is that, on each restart 
or iteration we can discard any angular momentum inconsistent, 
with a steady spin situation. Thus, we can avoid the problem 
of carrying along tumbling momentum which could considerably 
lengthen the computational time or cause the solution to 
pass right through a steady spin mode or recovery without 
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stabilizing. Of course to realize this advantage it is 
necessary to be able to differentiate between that part 
of the motion which constitutes tumbling and that part which 
is consistent with steady spinning. Herein lies one of the 
major advantages of writing the equations of motion in an 
inertial cylindrical reference frame rather than in a body- 
fixed moving coordinate system, for in the hybrid inertial 
system tumbling is readily recognizable as terms involving 
^ ^ or 0 . 

B. SIMPLIFICATION OF THE EQUATIONS OF MOTION 

With the above approach in mind the equations of motion 
can now be simplified sufficiently to meet the requirements 
of the solution scheme without the necessity of making any 
restrictive simplifying assumptions. Also since the 
intended iterative technique deals with the problem directly 
in dynamical units, we do not have to non-dimensionalize as 
would be necessary if one resorted to a somewhat more 
abstract mathematical method of varying the spin parameters 
to obtain convergence. As a result one can avoid obscuring 
the physical nature of the problem and is in a better 
position to utilize insight in monitoring, trouble shooting, 
and interpreting iteration trends. 

The simplification is started by recognizing that only 
steady or quasi-steady spin modes are sought. Thus, by 
setting angular orientation rates { S, ^ , and ^ ) equal to 
zero a considerable number of terms can at once be eliminated. 
The equations contained in Appendix A thus reduce to the 



following form: 
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Equation : 



= ~ (Cf^ 0(/i + CFy + Cp^ °^3s) ^ {(^^) ^ 



^ Equation : 



MR-Mt^R = - (Cp^ o(„ tCpy c^z, + Cp^ »<},) 

CR)^ ] ( 74 ) 



y Equation: 

m{zRURR)-^iy{ (ip*r)c(,3 i- e cos <j> tiy f(/^' 

— fc(Cpf^o(^3 ■^CpiyC{^ji-<^^o{j^J-R(C3yR/z 

■h Cpy 0 ^ 2 Z 
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( 76 ) 



^ Equation : 

I. ((H) ^ ^ cos<l>]cos4>-iyl(^^/)c(^^ 

' e$ifjf]sii>i<p -ix y^cos& s/f/<p^,3 

-Xyf^COS &CO$4 >o<23-<- S//J& C< 2 S 

= (c,,2<<^osf- eMy S/M<P) ^sc 

Equation : 

” ( ^Aiy^JS '*' ^/£ K?J ) '*'(^J ] 



0Equation ; 




■hly^ ^S//J e S///0 c/2 3 



( 78 ) 



Although R is a rate parameter which is not consistent 
with the steady spin assumption, it effects both the values 
of cK and ^ and also the magnitude of the relative wind, and 
therefore must be retained in the equations. 

Next, to put the equations into a useable form for the 
suggested iteration scheme we simply solve each equation 
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for the corresponding acceleration parameter in terms of the 
other variables. The equations then become: 

Z Equation: 



K Equation: 

^ ) ZM (00) 

}f Equation: 

jf= «';/ -t Kjj } *^^<13 

t ( C.(Cfjf^ CI^J + Cf^y <^23 ^M^33 )~^ (^^12 * ^ ‘^'22 

0 Equation : 

e ^(i,cos‘4>3lyS/2J^^]'[-(XrXrX'^^^)^‘^^‘2"J^cosi^ 
3-{l,$m^<p3-lyCOS^4i-Il)'i!^ SM&QOSe* Coif ( 82 ) 
-c„y Sii^(^] ] 
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Equation : 




sin<pcosj>-i^ 0O)S $+ k»y,p<i} 
1-Cny<>^l3+0^‘<S3] ^ [(/#+(Kf+(0*'J ] 



( 83 ) 



(j) Equation : 




^ ^ ^ s/fj d> COS (p -{fi-r)cose 



( 84 ) 



Before procee(3ing to develop the details of the proposed 
iterative scheme of solving these equations, some comments on 
the steady spin assumption should be made. First of all, 
truly oscillatory spin modes, as opposed to lightly damped 
steady spin modes that may be very slow in convergence 
should not, in general, result from the use of mean force 
and moment model test date, the usual form of such data. 

Thus if the oscillatory nature of a spin were attributable 
to periodic vortex shedding from the fuselage or stalled 
aerodynamic surfaces (von Karman Vortex Street phenomenon) , 
the use of mean aerodynamic forces and moments should 
effectively filter out such oscillation from the computer 
solution. If however the oscillations were strictly a 
result of the dynamics of the problem, such as a limit 
cycle type of oscillation, then the oscillatory nature of 
the spin should appear in the computations as a failure to 
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converge to steady spin parameter values. In this case, the 
complete equations would have to be used along with a very 

small time increment and the stability derivatives due to ^ 

* 

and ^ to obtain an accurate description of the motion in the 
form of a numerical time history. Although this would involve 
solving the full equations for the corresponding acceleration 
terms as was done in equations (79) through (84) and would 
also lengthen the computation time it would not be a diffi- 
cult extension and might well provide some valuable insight 
into the stability of such oscillatory modes. This form of 
the spin equations would be especially valuable in studying 
the effect of applying oscillatory recovery control inputs in 
resonance with the spin oscillations to obtain an "autopilot 
spin recovery" from an otherwise irrecoverable spin mode. 
Another use of the full equations would be to study the 
stability of spin solutions obtained by other numerical 
methods. However terms involving 0/ ^ / and 0 would re- 
introduce the effects of angular momentum into the computations 
and would probably cause excessive tumbling if the full 
equations were used to initially solve for spin modes. There- 
fore either the following computational scheme or some other 
numerical technique for solving non-linear differential 
equations such as the iterative procedure described on page 
270 of Ref. 7 should be employed to initially solve for spin 
modes . 

C. SUGGESTED COMPUTER SOLUTION SCHEME 

Although the equations are rather complex the suggested 

computer iteration scheme is fairly simple. Basically some 
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starting position or initial "guess" must be supplied in 
terms of altitude rate, spin radius , spin rate, orientation 
(©, ^ , and ^ ) , along with cont'rol positions ^ and ' 

Angle of attack and angle of sideslip are then computed from 
equations (65) and (66) and, in turn C><^ and ^ are used to 
obtain values for the various aerodynamic derivatives dis- 
cussed in Section IV B by interpolating from appropriate 
tables. The six aerodynamic force and moment coefficients, 
and^^/ are then computed as outlined in equations (71) 
and (72) and, with these values in hand, and X can be 
computed directly from equations (79) and (80) . The other 
accelerations , Q , y , and 0 , can be computed by a 
looping routine encompassing equations (81) through (84) . 

This looping will be necessary due to the coupling of the 
respective acceleration terms in each of these equations. 

Once consistent values of Y , Q > y > snd (p have been 
obtained from the looping iteration, incremental changes in 
each of the seven computation parameters can be obtained 
from the following equations; 

s' in ( 85 ) 






( 86 ) 






(87) 



^Radius rate, , is initially set equal to zero. 
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(88) 









( 89 ) 






( 90 ) 



" At’- 






( 91 ) 



If the full equations were to be used, obviously it 

^ ♦ 0 

would then be necessary to also compute e.ff ,4> . 

^ , and <p as follows: 



^rw+o 



( 93 ) 






( 94 ) 



‘A Af 






( 95 ) 



^(yiti) ~ P(n) 



( 96 ) 
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( 97 ) 



The computational sequence outlined above is then repeated 
using the revised values of the parameters involved and 
continued until some pre-designated convergence criteria is 
met, such as all accelerations being less than some small 
epsilon value or the rates of change of the various 
parameters are each in turn less than some small character- 
istic values as individually specified. For example, one 

2 * 

foot per second in /^ , 0.5 feet per second in , one degree 
per second in ^ , etc. 

The choice of an incremental time probably will have a 
very significant effect on the overall success of this 
iteration scheme. The use of a very small value of ^t.will 
of course provide the smoothest convergence but at a 
significant expense in computer time. A large At on the 
other hand can easily preclude convergence by drastically 
altering all parameters from iteration to iteration. To 
illustrate this point and also pin down what is meant by 
large and small At it is necessary to step back a bit and 
consider conceptually what is really being done by using an 
iterative solution scheme. As pointed out earlier, the 
equations are highly non-linear and extensively coupled, 
and consequently one has very little insight as to what 
linearizing approximations would be valid or, for that 
matter, over what range of values can the effects of an 
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individual parameter change be considered linear. This 
problem is however overcome in the suggested iteration 
scheme by allowing the dynamics of the problem to pick the 
magnitude of the incremental changes in each parameter. 

Thus if the At used is of such a magnitude as to constitute 
a small pertubation relative to the characteristic time of 
the problem, the time per revolution, we essentially estab- 
lish appropriately sized small pertubations in each of the 
spin parameters as well. Thus, for example, a spin rate of 
40° per second would yield a 9 second turn and possibly a 
good starting guess for an appropriately sized would be 

.09 seconds. The validity of this concept will not be 
rigorously or otherwise established here, other than to say 
that in general small pertubation techniques have enjoyed 
some success in other engineering applications and this 
appears to be a quite natural extension. 

The problem of picking appropriate starting guesses for 
^ , and (p unfortunately is not as easily 

handled as deciding on an initial/^t, . Although some rather 
broad guidelines can be stated such as steep nose down 
pitch, large radius and low spin rate for normal erect 
(upright) spins and shallow nose down pitch, small radius 
and high spin rate for flat spins, the actual searching for 
individual modes must proceed on a trial and error basis 
until an analysis of solutions suggests some better pro- 
cedures. Starting guesses that don't lead to spins will of 
course show up as recoveries to linear stalled flight atti- 
tudes or complete fly-away recoveries. 



43 



Certainly some insight can be used to pick as a 
negative (down) and of an order of magnitude consistent 
with stalled flight velocities (T.A.S.) at the altitude 
where the solution is sought. Likewise, a wings level 
attitude with the nose down and pointed at the spin axis, 
e.g. , 0 “ +30°, V ~ -90° and ~ +90°, would be a reason- 
able starting orientation for an erect right spin. Finally, 
0 

^ should obviously be chosen small so that it will not toss 
the airplane out of the ball park under the action of 
centrifugal force before the aerodynamics can come into 
play. One point should be made however concerning searches 
for flat spin modes with small radius values. If conven- 
tional force modeling is used, as discussed in Section IV B, 
then in addition to the general lack of validity of adding 
non-linear aerodynamic forces and moments, there is the 
compounding of errors due to the fact that free stream o( , 

^ and ^ vary from point to point over the fuselage and 
wings as a function of orientation when radius is reduced to 
the same order of magnitude as fuselage length or wing span. 
Thus, for example, if the wings are level and the nose is 
pointed at the spin axes in a fairly flat attitude the 
geometry of the problem will yield a significant difference 
in both sideslip velocity, and consequently from nose to 

tail for a small spin radius. Of course aerodynamic data 
for this type of flow condition cannot be obtained in a 
static wind tunnel setup. 
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D. AUXILIARY ANGULAR RELATIONSHIPS 



The particular set of "Euler Angles" used in the 
derivation of equations (79) through (84) is just one of 
several variations of such positioning angles which are 
generally classified as Euler Angles. The principal 
advantage gained by the use of such sets of positioning 
angles is that they describe a single fixed spatial 
orientation of a rigid body for given values of e, ijj . 
and (p , and yield a convenient set of equations for the 
direction cosines with the order of angular rotation 
implicit in the equations. However unlike the more con- 
ventional pitch, roll and yaw orientation angles which can 
be readily visualized once the order is specified, 
orientation in terms of Euler Angles is difficult to visualize 
without some sketching, which can be readily appreciated by 
refering to Fig. 2. Thus unless the orientation is close 
to quadrantal a computer output of orientation in Euler 
Angles may be difficult to interpret. 

To alleviate this situation, relationships between an 
ordered set of orientation angles (1) pitc±i, (2) roll, and 
(3) yaw, and the Euler angles shown in Fig. 2 were 
derived. The derivation was facilitated by the sketch of 
the spherical triangle relationships between the two 
angular sets shown in Fig. 4 and involves the use of 
Napier's analogies from Spherical Trigonometry, Ref. 5. 

The use of Napier's analogies avoids the quadrantal 
ambiguity associated with the solution of right spherical 



45 



triangles by other methods. Thus the equations for con- 
verting the ordered angular rotations Q, (!) and (Py to the 

P ^ T 



set of Euler angles used in this paper are: 






^^9* 4 Schematic of Auxiliary Angular Relationships 
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The equations for obtaining , 0^ and Euler angles 

are : 




Q 






Cos(1^^j 






CoS 



:io2) 



(j)-^^0-ZTan 



-I SlhJ 



.if 



|t4w'| 



5(W {iS^) 



i] 






(103) 



Equations (101) through (103) can be converted for computa- 
tions in radians by substituting for 90°. 
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VI. THE INVERSE PROBLEM 



As discussed in Section III B the conventional method of 
modeling the aerodynamic forces and moments on the spinning 

I 

airplane with static wind tunnel data lacks mathematical 
validity. Therefore it will be necessary to develop a 
better method of modeling the aerodynamics of the problem 
if the computational approach to predicting spin modes is to 
be successful. The work done with the rotary balance. Ref. 

4 has been one effort in this area. Such work, however, is 
difficult to pursue without some method of obtaining the 
actual forces and moments as a basis for comparison. 
Apparently in the past the validity of such aerodynamic 
data has been checked by the indirect method of using it to 
solve the partially linearized equation and then comparing 
computed spin modes with known spin characteristics. This 
approach may allow some estimate of the overall success of 
the combined result of the equations and data but doesn't 
establish the validity of either. A better approach would 
be to use the exact equations derived in Section V in such 
a manner as to allow comparison of the 6 computed aero- 
dynamic forces and moments with the corresponding forces 
and moments on the actual spinning airplane. Fortunately 
this can be done for the steady spin case and thus allow 
some progress to be made in modeling steady or mean values 
of forces and moments. 
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In order to use the equations of motion derived in 



Section V in this application it is first necessary to 
obtain the values of 2^ / K , ^ ^ , and ^ for the actual 



airplane or free spinning model. For the case of the model 
in the free-spinning wind tunnel discussed in Ref. 6, 
determination of the above spin parameters can obviously be 
obtained from tunnel speed, and photographic coverage, and 
thus will not be discussed further. For the actual airplane 
the means of getting these parameters is still fairly 

simple but in some cases less obvious, and consequently 

$ 

will be outlined here. First of all, can be determined 

by the pilot altitude/time observation or photo panel data, 
or externally by ground station precision radar or optical 
tracking. Likewise If can be determined either by the pilot 
or from ground station photographic coverage. Spin radius 
can be obtained by resolving triaxial eg accelerometer data 
along the principal axes and using the following relationships 




( 104 ) 




( 105 ) 



Airplane orientation in terms of Euler Angles ^ ^ , and 



^ can be determined from values of pitch, roll, and yaw 
rates and the triaxial accelerations as follows: 



? =- 



( 106 ) 
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if Si^e cos<p 



(107) 



r = ^ cos e 



Thus : 

^=7-.«Y— ^ ] 



(108) 



(109) 

(110) 



and from Table I : 
or solving for f7 .. 



(/;=cos 



rij 



yr*o%*»# ^^ >**>%*<» « ^*** ^ »**^^"*w K t u m u mnw 






( 111 ) 



T "'/ J^Slf>i 0COS^+^/CO$0COSO Si 

" %CoS(p-Sy 

Obviously there will be some loss of accuracy in the 
determination of spin radius and orientation if the spin is 
not exactly steady. However small fluctuations in the 
triaxial accelerations and pitch roll and yaw rates can be 
averaged out and should allow computation of representative 
values of , Q , ip , and for the near steady spin case. 

Next, equations (79) through (84) can be adapted to the 
problem at hand by setting R and all acceleration parameters 
equal to zero. The equations thus reduce to the following : 
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(112) 



2- Equation: 



~ '"'x “^>3 ^Fy <Fix} * 



^ Equation : 

^RN\ 

^ Equation: 



'^U ^F/ I ^ 3 1 



R ( C ^2Z Kji. ) 



(113) 



(114) 



0 Equation : 

_ 2 ( Jx 5w ^4> t ly COS ’’4) -If ) f @ COSO 

Csc{cefL)^t(i.tj ai5) 

— Ci^^CO$(p-CMyStt<}^ 

fp Equation: 



^/3 



“CT 



(116) 



(P Equation : 



z(ix-iY)y^siM^e stt3<i><x)s<l> 




(117) 



Equations (112) through (117) can also be cast into matrix 
format for direct computer computation as follows : 
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(118) 
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(119) 



where the fact that , as expressed by equation (116), 

has been used in equation (114) to show that Cp^=Q in the 

Yl 

steady spin. 

In addition to allowing the computation of the force and 
moment coefficients from steady spin parameters the results 
of this development show that and are not 

A, f 
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necessarily zero in the steady spin case, contrary to what 
one might have intuitively assumed. In fact, the steady 



spin situation where are all zero is a 

special case not likely to be encountered in airplanes of 
conventional design. To illustrate this point, equations 
(115) and (117) can be solved for the case where the moment 
coefficients are all zero with the result that: 






( 120 ) 



and 

it 



( 121 ) 



or simply : 

Thus it is seen that in a conventional airplane where it is 
unrealistic to have equation (122) satisfied, aerodynamic 
moments must be present to balance the gyroscopic moments 
and gyroscopically precess the inertial orientation of the 
airplane sufficiently to maintain a fixed orientation 
relative to the spin axis and radius vector. 
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VII. RESULTS AND DISCUSSION 



A. ANALYTICAL ASPECTS 

Although the primary intent of this effort was to 
develop another method of computing airplane spin modes 
from wind tunnel test data, the equations of motion derived 
in Section IV can also be used to advantage in a qualitative 
analytical approach to the problem of determining spin 
characteristics. In addition to possibly suggesting new 
parametric indices of spin characteristics, these equations 
should also allow a number of experimentally observed effects 
of parameter variations to be put on a firmer analytical 
foundation, such as the effect of the sign of (I -I ) on 

^ X y 

spin recovery [Ref, 8] . The fact that the inertia character- 
istics appear so prominently in equations (79) through (84) 
confirms test results which have suggested that the gyro- 
scopic aspects of the problem are as important as the 
aerodynamics. Another result of this analysis of the spin 
problem was to emphasize the non-linearity of the aero- 
dynamics associated with the spin. Thus computational 
techniques which use superposition in modeling the aero- 
dynamics or employ partially linearized equations of motion 
should be expected to yield rather poor results where high 
spin rates or small spin radii are involved. These are 
probably the main reasons why the conventional approach has 
not really been successful in treating spins in general. 
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and has been particularly unproductive in predicting flat 
spin modes. 

B. THE COMPUTER PROGRAM 

Although the suggested computer program outlined in 
Section V C uses conventional aerodynamic modeling which is 
not entirely representative in the spin, it is felt that 
development of a computer program capable of solving the 
equations of motion is the next essential"” step in this 
approach to the spin problem. Once the computer program 
is working, improvements in the modeling of the forces and 
mements can be incorporated as they are devised. Thus the 
computational and aerodynamic parts of the problem can be 
pursued, in the interim, as separate problems. It is with 
this in mind that the author has chosen to discuss a few of 
the aspects of the computer solution problem rather than 
address the aerodynamics of the problem. 

The iterative scheme outlined in Section V B may exhibit 
a tendency toward a lightly damped oscillatory form of con- 
vergence since it is essentially following the dynamical 
motion of the airplane, which is characteristically lightly 
damped in a spin. The oscillatory nature of the convergence 
is compounded somewhat because we cannot be sure how far a 
starting guess may be from a force and moment equilibrium. 

To improve this situation there are a number of ways to 
provide some synthetic damping to the iteration to speed 
convergence. Some of these ways are: 
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a. If excursions in radius prove to be a problem as 
a result of momentum in radius rate (^) , the motion can 
be damped as follows; 

V/) = i % * 

where ^ is some factor between zero and one. 

b. If excursions in altitude rate become a problem they 
can be handled as follows; 

'^ 27 ) 

Another problem that may arise is that convergence in 
orientation may be deadbeat, depending on the size of At. . 
This is due to the fact that we have eliminated orientation 
rates to reduce any tumbling tendencies. Thus it may be 
necessary to augment the rate of change in orientation per 
iteration to speed convergence. This can be done in the 
following manner: 



' ^<») ^ ^ ^ 
■%*>) = ^ ( ^ 4 ^) 



(125) 



(126) 



(127) 



Where b is some factor greater than one. 

Finally, since our starting guess in ^ may be out of the 
ball park and cause large changes in radius prior to reaching 
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a magnitude consistent with the aerodynamics of the problem, 
it may be necessary to augment the rate of change of this 
characteristically unresponsive parameter. This may be 
accomplished in the following way: 



where C is some factor greater than one. 

C. FOLLOW-ON STUDIES 

As mentioned earlier, the equations of motion as derived 
in this paper should make it possible to provide an analyti- 
cal basis for many of the heretofore experimentally observed 
characteristics of airplane spinning motion. Such a cor- 
relation study would be relatively inexpensive to pursue and 
would probably provide a wealth of additional insight into 
the spin problem as a whole and possibly be beneficial to 
the more immediate goals of computer program development and 
force-moment modeling. It might also be instructive to 
program the full equations of motion, as contained in the 
appendix, on an analog computer having a digital computer 
interface capability for the aerodynamic data. Some other 
aspects of the spin problem which could be pursued once the 
computer programming and force modeling problems are solved 
are : 

(1) Determine the spin modes of currently operated 
tactical military airplanes and investigate the effective- 
ness of recommended spin recovery techniques. 




(128) 
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(2) Determine what characteristic aerodynamic parameters 
can be most economically fixed in existing airplanes with 
irrecoverable spins and investigate the magnitudes of such 
changes that would be required to make recovery possible. 

An example would be providing rudder stop over-rides or 
additional deflection for spin recoveries. 

(3) Study the feasibility of providing strap-on devices 
such as a deployable fabric ventral fins or wing mounted 
retro-rocket pods for use on tactics training flights. The 
size of the device or necessary thrust requirements could 

be fairly easily determined with the assistance of a working 
computer program, or, in the case of the retro-rocket, could 
also be determined by the "Inverse Problem" Technique out- 
lined in Section VI. 

(4) Determine "design envelopes" that prescribe limits 
in terms of common aerodynamic derivatives in order to 
assist the designer in avoiding troublesome combinations 
wherever feasible. Such work would also provide the basis 
for improvements in the current Military Specification 
MIL F-8785 (ASG) . 

(5) Study the feasibility of utilizing an autopilot 
to apply oscillatory control inputs in resonance with the 
dynamics of the spin to break an otherwise irrecoverable 
mode . 

These are only a few applications that might be under- 
taken to alleviate the spin problems inherent in the 
designs of todays high performance tactical military 
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airplanes. However the satisfactory accomplishment of just 
a few is all that is really needed. 
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VIII. CONCLUSIONS AND RECOMMENDATIONS 



The approach to the spin problem outlined in this paper 
should provide a means of obtaining some additional insight 
into the problem of aircraft spin. Once a successful 
computer program has been written and a more representative 
means of modeling aerodynamic forces and moments has been 
devised it should be possible to solve the equations of 
motion and obtain some idea of the spin modes that a given 
airplane design may exhibit. The results of Section VI 
indicate that gyroscopic moments play an important role in 
the spin problem. Thus it may also be worthwhile to treat 
the spinning airplane as a slow turning free gyroscope, of 
relatively high inertia, being acted upon by triaxial aero- 
dynamic forces and moments. Such a formulation of the spin 
would allow the application of the broad scope of knowl- 
edge of gyroscopic motion from the fields of physics and 
mechanics to the airplane spin problem. In this manner it 
may be possible to obtain some additional understanding of 
the stability of various spin modes and possibly avoid some 
undesirable airplane spin characteristics without having 
detailed knowledge of the aerodynamics involved. 
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APPENDIX 



SPIN EQUATIONS 

The following equations constitute the full equations of 
motion for a spinning airplane as written in a cylindrical 
coordinate system. 
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